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Abstract 



We use both old and new theoretical developments in QCD dispersion rela- 
tion constraints on the scalar form factor in the decay K — > irlvi to obtain 
constraints on the strange quark mass. The perturbative QCD side of the cal- 
culation incorporates up to four-loop corrections, while the hadronic side uses 
a recently developed parameterization constructed explicitly to satisfy the 
dispersive constraints. Using chiral perturbation theory (xPT) as a model for 
soon-to-be measured data, we find a series of lower bounds on m s increasing 
with the accuracy to which one believes %PT to represent the full QCD result. 
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I. INTRODUCTION 



The level of interest in using dispersion relations to study the analytic properties of 
Green functions has ebbed and flowed for decades. Recently, attention has increased due to 
studies of bounds obtained on form factors from heavy hadron semileptonic decays, which 
can be used in concert with Heavy Quark Effective Theory to isolate the CKM elements 
\Vcb\ and \ Vub\ (see [|l] for a compilation of references). Such studies were initially motivated 
by the 1981 paper of Bourrely, Machet, and de Rafael 0, which first applied the basic 
method of equating the dispersion integral over hadronic form factors to its perturbative 
QCD evaluation in the deep Euclidean region, in that paper for the case of K e3 decays. 
Here we revisit decays, armed with new improvements in technique and the promise of 
vastly more data to be produced at DA$NE [Q. Because the dispersive bound uses QCD 
to constrain the shape of each hadronic form factor (as a function of momentum transfer), 
one can use the experimentally measured form factors to obtain limits on QCD parameters 
such as the quark masses; however, until this data is available, one can use precisely the 
same method to take the form factor from a given model or calculation as being correct, 
and discover limits on values of the QCD parameters for which the shape of the form factor 
is consistent with the dispersive bounds.Q 

In particular, we use the one-loop corrected chiral perturbation theory (xPT) result of 
Gasser and Leutwyler |4j] for the scalar form factor in decays. The scalar form factor 
is especially interesting because it contributes to the scalar current correlator, which is 
proportional to (m s — m u ) 2 (i.e., vanishing for flavor-conserving processes), and is thus more 
sensitive to quark mass values. Moreover, it satisfies low-energy theorems in to a 

high degree of precision. Our results may also be used for the vector form factor to obtain 
interesting information, but our thrust in this paper is to show how experiment will be able 
to teach us about limits on m s , and specifically, that the consistency of and QCD 
places substantial lower limits on the strange quark mass. 

This paper is organized as follows: In Sec. II we summarize the most important parts 
of the dispersion relation approach with references to papers explaining the fine points in 
more detail, as well as describe our x?T inputs. Section III presents the QCD result for the 
scalar current correlator to four loops, including leading nonperturbative effects. In Sec. IV 
we present improved versions of the analytic functions used to parameterize form factors 
consistent with QCD and describe how to judge the quality of fits to this form. In Sec. V 
we present our results and conclude. 

II. USING THE DISPERSION RELATIONS AND %PT 

By means of Cauchy's theorem, the dispersion integral simply relates an integral of a 
Green function (in our case, a current two-point correlator) expressed over a kinematic region 



Alternately, one could insist that particular QCD inputs are correct and use the bounds to obtain 
limitations on the consistency of the model with QCD. Since we are using the model as a proxy 
for real data, we do not follow this approach here. 
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of momentum transfer q 2 where hadronic quantities conveniently describe the physics, to the 
Green function itself evaluated at a point deep in the Euclidean region, where perturbative 
QCD is reliable. In the present case, the current is = s'y^u, whose divergence d ■ V = 
i(m s — m u )(su) is just the scalar current multiplied by an explicit factor of quark masses; it 
is this feature that leads us to interesting bounds on m s . The scalar correlator, related to 
the vector correlator by means of a Ward identity, is given by 



ijj(Q 2 ) = i J d^xe iqx (0\Td fl V^(x)d l/ V u (Oy\0), (2.1) 
where Q 2 = —q 2 . The unsubtracted dispersion relation reads 



The n-times subtracted form of ( |2.2p is readily obtained by taking Q 2 derivatives on both 
sides, and in the remainder of this section we take ip(Q 2 ) to mean the perturbatively com- 
puted two-point function with any fixed number of subtractions. 

The l.h.s. of (|2.2j ) is evaluated at a point Q 2 far from the resonant region of V 1 , which for 
s — > u transitions is satisfied when Q 2 ^> Aq CD . In practice, we choose Q 2 = 4 GeV 2 , which 
is popular in lattice simulations. The r.h.s. of ( |2.2| ) is expressed as a sum of integrals over 
hadronic quantities. It is a manifestly nonnegative quantity, and so neglecting any subset of 
contributions produces a strict inequality, in that any subset of hadronic form factors in any 
kinematic region connected by crossing symmetry is limited by an expression computable 
in QCD. In particular, integrals of the form factors for K — > tt in the crossed kinematic 
region of (vacuum — > Ktt) are bounded, leading to restrictions on the shape of form factors 
permitted by QCD. The tightness of these bounds is of course regulated by numerical inputs 
to the QCD calculation, and particularly in this case by m s . It is precisely this sensitivity 
that is useful to us: Given the shape of a K — > tt form factor obtained from data or a model, 
one may ask which values of m s permit this shape to be consistent with QCD. 

The derivation of the bounds on the shapes of semileptonic decay form factors and its 
expression in terms of a well-defined parameterization draws upon techniques developed in 
Refs. f2|, and |5|-|7j, and the approach is explained in detail in ||. Here we present expressions 
only to establish notation for later use. Beginning with the kinematic points 

t± = (m K ±m n ) 2 , (2.3) 

it is convenient in such decays to use the kinematic variable 



where t s < t+ is a kinematic scale fixed for convenience, as described below. For each 
form factor F(t) there is a computable function <pF(t;t s ;Q 2 ) obtained from ( |2.2|) , and in 
terms of these quantities all functional forms for the form factors allowed by QCD may be 
parameterized by 



F(t) = -7 \, ; ' E*nz(t;t s r , (2.5) 
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where the coefficients unknown parameters obeying 

oo 

E M 2 < 1 • ( 2 - 6 ) 

n=0 

This result employs analyticity of the two-point function away from hadronic thresholds, 
crossing symmetry between the matrix elements for Ktt production and K — > tt, and knowl- 
edge of QCD at Q 2 . To lowest order in Gp this form is exact, for there are no resonant 
poles, multiparticle continuum states, or anomalous thresholds in the scalar channel below 
the threshold t = t + of (vacuum — > Kir). To be precise, one must also distinguish the two 
isospin channels (vacuum — > K°t[~) and (vacuum — > K~ir ). Apart from a small explicit 
isospin breaking in the form factors, whose inclusion is described below, the two channels 
have slightly different values of t±. Since nonanalytic q 2 behavior in the correlator, which 
arises from physical and anomalous particle thresholds, is what determines the form of the 
parameterization fl2.5| ), we take t + to be the smaller of the two possibilities, (%+ + m w o) 2 , 
so that no threshold is ever crossed in the region where ( |2.5|) is used. 
The two form factors relevant to decays are defined by 

(-K(p')\V»\K(p)) = 2f + (q 2 ) (V - ^V) + d(q 2 )^ (2.7) 

where g M = (p — p') M . This decomposition separates the vector / + and scalar d form factors, 
and consequently the bounds studied in this work apply to d. In the differential width, 
\d(q 2 )\ 2 appears multiplied by a helicity suppression factor of mj. 

In lieu of very precise experimental information, we use the results [|J of one-loop chiral 
perturbation theory for the scalar form factor as "data". At two kinematic points in par- 
ticular, q 2 = and q 2 = {m 2 K — m 2 ), predictions of amazing precision have been made. At 
q 2 = 0, the authors of || find 

ff n ~(0) = 0.977, /+ +7r °(0) = 0.998. (2.8) 

By the Ademollo-Gatto theorem the symmetry relation /+(0) = 1 is violated by terms 
of order m 2 , and we see these percent-level deviations in the predictions ( |2.8|) . Uncertainties 
on (|2.8|) are even smaller, perhaps a few parts per 10 3 . The relation to the scalar form factor 
at this point is 

d(0) = (m 2 K -ml)U(0). (2.9) 
At q 2 = {m 2 K — m 2 ), the value of d(q 2 ) is given by the Callan-Treiman relation 0, 

d{m\ - ml) = (m 2 K - m 2 n ) + A CT ^j , (2.10) 

where fic/fiv — 1-22 ± 0.01, and Act measures deviations from the exact Callan-Treiman 
limit. As calculated in |3J], Act = —3.5 ■ 10~ 3 . Since no strange + resonances appear 
until the relatively high mass fQ(1430) and Kq(1950), the form factor in the decay region 
should be very smooth (and indeed, almost linear) to a very good approximation. It is 
therefore very reasonable to assume that the scalar form factor is given in the intermediate 
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region, i.e., between q 2 = and q 2 = {m? K — m 2 ), by chiral perturbation theory to the 
same accuracy as at the endpoints. It is, of course, difficult to estimate this error, although 
some indication may be obtained from a forthcoming two-loop xPT calculation [ID], despite 



the possible appearance there of currently undetermined renormalization constants. At any 
rate, since the one-loop corrections to the low-energy theorems are already quite small, the 
whole framework of would crumble if new experiments showed deviations from the 
predictions by more than a few percent. 



III. THE SCALAR CORRELATOR IN QCD 



Instead of the divergent scalar correlator ip(Q 2 ), we consider its second derivative ip"(Q 2 ), 
which is free of renormalization point dependence and therefore satisfies a homogeneous 
renormalization group (RG) equation. ifj"(Q 2 ) is calculated in perturbative QCD for large 
Q 2 using an expansion in powers of the quark masses and the operator product expansion. 
The leading term has been calculated to four loops ]nj] , and the O {m 2 / Q 2 ) correction to three 
loops |T2]]. The 0{m i s /Q i ) contribution is best considered together with the nonperturbative 
condensate terms. Collecting these results, 



6(m s - m u ) 
(4vr) 2 Q 2 
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7T 
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where the RG-invariant condensate combinations I s and Iq are given by 
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(3.2) 



(3.3) 



where we have set rif = 3 and omitted logarithms that vanish when taking fi = Q = \fQ 1 . 
Note that the peculiar 7i/a s terms cancel. The full result depends logarithmically on the 
renormalization point \i and on the parameters of the theory, like a s , m s , and condensates, 
which are renormalized at /i; however, as has been advocated in [0 , we implement the RG 
improvement for the case of the scalar correlator in the following way: ip"{Q 2 ) is evaluated 
at fi — Q, and the parameters a s and m s are extrapolated from a chosen reference point 
(in our case A^) to /i = Q using the four-loop beta functions (compiled in |L3|]). The 
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condensates are so poorly known and their effect at the chosen scale Q 2 = 4 GeV 2 so small 
that their \i dependence may be ignored. 

The numerical values chosen for the QCD inputs are Q 2 = 4 GeV 2 , A^~ = 380 ± 60, 
m u = m s /25, (m s uu) = (m s ss) = -f 2 K m 2 K = -0.031 GeV 4 , and (a s G 2 /n) = 0.02-0.06 
GeV 4 , although the sensitivity of the analysis to these nonperturbative parameters is small. 
The scale Q is arbitrary subject to the constraints that if it is too small, perturbative QCD 
is unreliable, while if it is too large, the dispersive bounds thus obtained are weak. 



IV. PARAMETERIZATION AND QUALITY OF FIT 



In order to use the parameterization of Eq. (|2.5|) , one requires expressions for the function 
for each form factor. In notation designed to be similar to that of [0], we define 



t- 



t. 



02 



lt+ + Q 2 
ti tf. 



(4.1) 



Note that, in our notation, one must set the parameter t s = —Q 2 to obtain the limit of 
Ref. P, in which case 02 — 1- The for each form factor defined in Eq. (|2.7|) is given by 




(3i + 
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01 + 
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l+Z 
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(4.2) 
(4.3) 



These expressions for differ from those in 0, because the dispersion integrals are formu- 
lated using different linear combinations of the two polarization tensor component functions 
than used in the other work, leading to a different pattern of subtractions. The factors rj 
represent Clebsch-Gordan coefficients based on isospin symmetry; that is, both charged and 
neutral K processes contribute to the dispersion relation, so we exploit the near equality of 
their contributions. For the decay K L — > 7i ± £ T i / £ (whose width equals, by CPT, either that 
for K° or K° semileptonic decays), rj 2 = 3/2, while for K + —>■ 7i°£ + b>e, rj 2 = 3. If isospin 
breaking between charged and neutral K processes is significant, one can incorporate this 
difference into rj 2 , or even include a factor with q 2 dependence (written in terms of z) if the 
pattern of isospin breaking is known. In practice, we obtain a conservative correction to 
77^-+ = 3 by supposing that the form factor ratio represented by ( ^.8|) at q 2 = persists for 
all values of q 2 in the production region, to obtain an effective 77^-+ 



2.92. 



In the original work [[|, the bounds are expressed as the determinant of an n x n 
semipositive-definite matrix, where the form factor is assumed known at (n — 2) points. 
This produces an envelope of allowed form factors resembling a chain of sausages, since the 
form factor is required to pass through each of the (n — 2) points. Although presented explic- 
itly in for only 1 or 2 fixed points, the determinant method can of course be generalized to 
an arbitrary number of points, with a corresponding increase in the complexity of algebraic 
expressions appearing (see, e.g., [0). However, once the envelope of points allowed by the 
dispersion relation and chosen fixed points is established, it is not true that any curve lying 
within this envelope still satisfies the dispersive bound. In contrast, the parameterization of 
Eq. (|2.5|) subject to (|2.6|) always satisfies the determinant bound of arbitrarily high degree. 
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The central features that make the parameterization Q2.5|) useful are the bound ( |2.6| ) 
on the coefficients a n and the smallness of the kinematic variable z over the full range for 
allowed semileptonic decay. It follows that one may express the form factor over the full 
range using only the first N + 1 parameters {oq, a%, . . . , a^}, with the remaining infinite set 
bounded in magnitude and forming a theoretical truncation error M, 6^: 



N \z ]N+1 



EH 2 ' ; = , (4.4) 



which means that the form factor fit to these parameters using ( |2.5| ) has a theoretical un- 
certainty no larger than 5n- 

The kinematic parameter t s is used to minimize the size of this already small truncation 
error Equation (4.4) makes it clear that this minimization occurs when z = lies 



within the kinematic range chosen for the fit, t G [t m m,t max ]] from ( |2.4| ) one sees that this 
can occur only if t s also assumes some value in this range. Therefore, the truncation error 
is minimized by plotting ( |4.4j ) as a function of t, t s G [t m in, t mauX ], and finding that value of t s 
for which the maximum value over all t G [t m in, ^max] is smallest. 

Now suppose that the "data", in our case the xPT expression for the form factor, has 
an "experimental" uncertainty A. In order to state that %PT "data" agrees with QCD to 
within A, it must be possible to expand the QCD fit to an order N such that Y<n=o \ a n\ 2 < 1 
and 5n < A. In other words, in order to be certain that the exact all- orders QCD form 
factor -F exact (which is not accessible to us) lies within an uncertainty A of the data F x p^, we 
require that the deviation 6~n of the truncated QCD form factor F trunc from F exact is smaller 
thanP] the experimental uncertainty A. This is the central argument of our reasoning. 

The requirement that both 5n < A and Q2.6| ) are satisfied is the key to obtaining bounds 
on parameters in ip"(Q 2 ), in particular m s . A curve such as F x p^, which does not a priori 
satisfy fl2.6|), when expanded as a power series in z eventually produces a value of a at for 
some iV so large that ( |2.6| ) is violated. This does not necessarily mean that F xPT violates 
QCD, because it is only required to equal the -F exact within an uncertainty of A. However, 
the value of N must be large enough that 5n < A, or else F trunc expanded to n = N is 
not necessarily a good enough approximation to use in place of -F exact ; if it is not possible 
to carry out this fit and maintain Eq. ( |2.6| ), then the chosen inputs are not consistent with 
QCD bounds. When the fit is successful, one concludes that the true QCD form factor is 
actually given by the fit of the form factor to the parameterization up to order N, 
and the higher-order terms cannot be probed with %PT since they lie within our stated 
uncertainty A. 

The particular nature of the fit is irrelevant to us, since one requires additional physical 
input to distinguish two curves that otherwise satisfy all the requisite conditions. For exam- 
ple, we use (|2.5|) and choose to Taylor expand the xPT form factor in z as an expression for 

F(z)(pp(z)/ \/^j"(Q 2 ), to obtain a n values until ( |2.6| ) is violated and throw away all higher 
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If one prefers to combine the theoretical and experimental uncertainties, then the requirement 



becomes |F tmnc - F xPT | < J A 2 - 8%. 
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a n . Alternately, one may choose a highest-order N and perform a \ 2 fit to {do, ai, . . . , a^v} 
subject to the constraint (2.6); this would be the preferred method of analysis with binned 
data. 

Now consider the dependence of the bounds on the strange quark mass, which enters 



approximately linearly in yip"(Q 2 ) as seen in ( p.l|) . If one fixes all other quantities and 
decreases m s , it is clear from (|2.5| ) that one obtains the same form factor by increasing each 
a n by the same factor. However, eventually the large enough to saturate (p.6|) , and 

thus one finds a minimal allowed value for m s . There is no corresponding maximal value, 
for one can certainly make each a n as small as desired and still satisfy Q2.6|) . Furthermore, 
as we increase the precision A to which we believe the xPT form factor holds, Sn must 
decrease, so that we must go to a higher order N in the parameterization expansion. It 
becomes increasingly more difficult to avoid large values of a n unless m s is increased, so the 
lower bound on m s becomes larger as A decreases. 



V. RESULTS AND CONCLUSIONS 

We present lower bounds on m s = m^ s (l GeV) from the form factor F(t) = d{t)/{m 2 K — 
m 2 ) of K + — > 7r°£ + z/£, since its bound (due to the isospin factor rj) is roughly y/2 tighter than 
that from decay, as described above. The natural tendency for the truncation errors 8^ 
is to decrease an order of magnitude with each unit of N, since 2 max ~ — z min ^0.1 for 
optimal values of t s over the range t min = 0, t max = (m 2 K — m^); however, this is somewhat 
counteracted by the fact that 5n oc y/W' ', which increases when the smallest allowed value 
of m s increases with N. Moreover, 5n can be much smaller than the numbers we present, 
since we do not include the effects of the \Jl — J2n \ a n\ 2 term in ((4.4p; near the saturation 
of ( |2.6|) by the first few terms, this is a large suppression. 

Our numerical results are summarized in Table I. The conclusion we draw is that, if the 
xPT form factor calculation is believed only to a precision of 2-5%, then we conclude only 
that m s > 40 MeV. At a precision of ±-1%, m s > 90 MeV, and at 1/20%, m s > 140 MeV. 
The same sort of analysis will be possible with data, although the details of the fit will differ 
somewhat. 

We conclude by comparing briefly with other recent approaches |T5|-|r7] that bound 
m s using dispersion relations. Typically, what is done is to saturate the hadronic side as 
much as possible with phenomenological input from expected pole contributions, and/or a 
continuum contribution in the deep Minkowski region of the hadronic integral modeled by 
the perturbative QCD result. Often, multiple constraints are obtained by taking moments 
of both sides. We view this work as complementary to those approaches. On one hand, it 
is minimal in the sense that only Kir, which is presumably but a small portion of the total 
hadronic result, is used in the bounding inequalities; of course, modelers may obtain stronger 
bounds by including Kq poles or continuum contributions at the cost of introducing model- 
dependent inputs. Although we considered only one moment in our calculation, certainly 
additional moments also give constraints, although then Q 2 must be adjusted to larger 
values in order to make sure that ip^(Q 2 ) is calculable perturbatively. On the other hand, 
it includes data from semileptonic form factors measured over their entire kinematic range, 
and so is expected to provide a substantial amount of additional input to such constraints. 
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We expect that combined program of calculations from all of these approaches will deliver 
rather strong constraints on m s . 



N 


{a n } 


t s (GeV z ) 


ip" ■ 10 5 


max 


AF 

L -* ± max 


m s (MeV) 


1 


a = +0.910 
ai = +0.414 


0.144 


1.16 


1.4 ■ 10~ 2 


2.8 • 10~2 


> 41 


2 


a = +0.418 
at = +0.184 
a 2 = -0.890 


0.142 


5.54 


3.1 • 10~ 3 


3.4 • 10- 3 


> 90 


3 


a = +0.273 
a x = +0.118 
a 2 = -0.584 
a 3 = -0.756 


0.140 


13.01 


5.0 ■ 10~ 4 


1.5 ■ 10" 4 


> 139 



Table 1. Bounds on m s = (1 GeV), fit parameters, and truncation errors based on the 



- K 



X PT|- 



saturation ( |2.6|) of the parameterization (\2.3j) . AF is an abbreviation for \F tl 
The numerical value of 5 n, max here neglects the ^/l — J2n \ a n\ 2 term in (\j.4j), which can be 
a large suppression when (as here) ( |2.6| ) is nearly saturated by the first N terms. 
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